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ABSTRACT 

Despite the ubiquity of M dwarfs and their growing importance to studies of exoplanets, Galactic 
evolution, and stellar structure, methods for precisely measuring their fundamental stellar properties 
remain elusive. Existing techniques for measuring M dwarf luminosity, mass, radius, or composition 
are calibrated over a limited range of stellar parameters or require expensive observations. We find 
a strong correlation between the ATs-band luminosity (Mr), the observed strength of the /-band 
sodium doublet absorption feature, and [Fe/H] in M dwarfs without strong Ha emission. We show 
that the strength of this feature, coupled with [Fe/H] and spectral type, can be used to derive M 
dwarf Mk and radius without requiring parallax. Additionally, we find promising evidence that the 
strengths of the /-band sodium doublet and the nearby /-band calcium triplet may jointly indicate a- 
element enrichment. The use of these /-band features requires only moderate-resolution near-infrared 
spectroscopy to provide valuable information about the potential habitability of exoplanets around 
M dwarfs, and surface gravity and distance for M dwarfs throughout the Galaxy. This technique has 
immediate applicability for both target selection and candidate planet host system characterization 
for exoplanet missions such as TESS and K2. 

Subject headings: stars: low-mass—stars: fundamental parameters—techniques: spectroscopic—stars: 

abundances—stars: activity— planets and satellites: fundamental parameters 


1. INTRODUCTION 

M dwarfs account for the majority of nearby stars 
(e.g. lRcid et al.ll2002i) . but due to their relative faintness 
and complex molecular spectra, precise (< 10%) mea¬ 
surements of their stellar parameters remain challenging. 
Measurements of the stellar mass, radius, luminosity, and 
composition of M dwarfs are important for studies of 
stellar structure ( Torres et al.ll2009h . Galactic evolution 
(jBochanski et all 1201. 'ill , and particularly for studies of 
exoplanets, since indirectly derived planetary character¬ 
istics and histories depend on accurate knowledge of the 
host star. M dwarfs are attractive planet search targets 
because their smaller sizes and lower luminosities result 
in larger radial velocity (RV) and transit depth signa¬ 
tures compared to Sun-like st ars. T hey appear to host 
n ume rous planetary systems ([Dressing fc Charbonneaul 
120151) and are prioritiz ed in planned exoplanet surveys 
le-g. iRicker et al1l2015f) . However, the characterization 
of the result ing planets and whether they are in the hab¬ 
itable zone (iKastind I1993T) could be clouded by lack of 
efficient and precise stellar characterization. 

The complexity of M dwarf spectra encumbers spectro¬ 
scopic modeling and the measurement of stellar param¬ 
eters. Numerous techniques have been developed to cir¬ 
cumvent this difficulty, by relating M dwarf parameters 
to easily-measured photometry or moderate-resolution 
spectroscopic features. These relations are calibrated 
using stars for which the parameters are known a-priori, 
such as M dwarfs with interferometrically-measured radii 
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(IBovaiian et all [20121 : IMann et al.l 120 1 . 311 : INewton et all 
I2014bli or M dwarf companions to high er-mass stars 
with well-measured c o mpos i tions (e.g. iBonfils et al.l 
20051: IJohnson &; AppsI 120091: iRoias-Avala et al.l 120121 : 
Mann et al.ll2013aHNewton et al .1120Mall . 

Determinations of M dwarf radius, luminosity, and 
mass can be made through empirical relations between 
these parameters and M dwarf effective temperature 
(T e ff), which can be estimated from moder a te-resolution 
spectr al features ([Bovaiian et all 120121 IMann et all 
I2013bli . However, this technique is limited by systematic 
errors when using stellar models t o predict M dwa rf pa¬ 
rameters from T e g and [Fe/H] (e.g. lBovaiian et aLll2012lh 
by the paucity of interferometric calibrator stars, and 
by a large range of M dwarf paramete rs at a given es¬ 
timated T e ff fe.g. INewton et al1l2014bl h Combining im¬ 
proved T e ff measures and [Fe/H] enables better M dwarf 
par ameters, but requir es broadband flux-calibrated spec¬ 
tra (IMann et al.ll201 5f). 

One can also estimate M dwarf parameters from K s - 
band luminosity (Mk), which is well-corre lated with stel¬ 
lar m ass and insensitive to mctallicity (Delfosse et al. 
1200011 . This allows mass measurements of nearby 
M dwarfs with parallaxes, but parallaxes for many 
nearby M dwarfs may not be available for several years 
dPerrvman et al.l 1200111 . For M dwarfs without parallax 
measurements, some success has been achi eved i n deter¬ 
minin g Mk using high-resolution spectra dPineda et al.l 
I2013IL but such spectra are observationally expensive. 
Methods that do not require parallax, high-resolution 
spectra, or absolute flux measurements of the target 
stars, and which can be applied immediately are highly 
desirable. 

We have developed a simple, powerful technique, based 
on the strengths of the /-band neutral sodium doublet 
(~ 820 nm, hereafter Na i) and singly-ionized calcium 
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triplet (~ 860 nm, hereafter Ca n), to enable deter¬ 
mination of M dwarf Mk and radius, and convey in¬ 
formation about M dwarf surface gravity and possibly 
a-enrichment. We present here an observational study 
that shows the potential applicability of these /-band 
features in conjunction with measurements of M dwarf 
T e fj and [Fe/H], all obtained from moderate-resolution 
near-infrared spectra. 


2. OBSERVATIONS AND METHODS 

The 342 stars used her e were observed in our IRTF- 
SpeX fRavner et aLil2003j ) M dwarf survey. iTerrien et alJ 
(|2012t) detailed the observations and data reduction. 
We also inc l ude spe ctra of nine nearby M dwarfs from 
IMann et alJ (|2013bfl . All spectra have R ~ 2000 and 
span 0.8 — 2.4/mi. We measured [Fe/H] for each star us¬ 
ing the weighted (by signal-to-noise, S/N) average of the 
J, //, K s NIR spectral cali brations offMann et alT(l2013at) 
(for M5 and earlier) and IMann et a0 ~( 2014fl (later than 
M5), which have <0.1 dex precision. We also measured 
the T e fj-sensitive H 2 O-K 2 index and associated NIR M 
spectr al ty pe (NIR SpT) explore d in iRoias-Avala et all 
(|2012t) and lNewton et all (|2014a[) . 

For each target, we compiled parallax (from 
the RFCON^ n dat abase van Altena et al.l 
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(iZacharias et al. 1 201 .'ll) and RV (|Chubak et ahl l2012t 
iNewton et al.l l2014afl where availabl e. We used Ha 
equivalent width measurements (e .g. iGizis et ah] 120021 : 
iRiaz et aJll2006l iGaidos et al.ll2014h to constrain activity 
levels for our M dwarfs. 

Our technique relies on measurements of the pseudo¬ 
equivalent widths (EWs) of the Na 1 and Ca 11 fea¬ 
tures. The N a 1 doublet is a known gravity indicator 
for M dwarfs jSchiavon et al.l 119971 : IMartin et al.l 1201 (1 
ISchlieder et al .1120 1 2li and both the Na 1 doublet and the 
Ca 11 triplet ha ve shown possi bl e sensitivit y to Ha-based 
activity levels ([Kafka fc Honevcuttl 120061 ). The defini¬ 
tions of the features and continuum regions are shown 
in Figure [T] For the Na 1 doublet and each component 
of the Ca 11 triplet, we fit a line to the two neighboring 
continuum regions (Figure [T|), and calculated the EW of 
the feature relative to this line. Our final EWc a value is 
the sum of the three components. We define positive EW 
as absorption for the Ca 11 and Na 1 features. The mea¬ 
sured and literature parameters for each star are listed 
in Table [T| For six M d warfs observed in our program 
and IMann et al.l (|2013bl ). we find no offset for the Na 1 
feature, but that our measurements of EWc a are system¬ 
atically 0.96 A higher, possibly due to our use of SpeX 
(R ~ 2000) in the /-ban d instead of SNIFS (lLantz et al.l 


12004 R ~ 1000) used in IMann et ah (12013b ). W e app l 
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this correction for the 9 stars from IMann et al.l (l2013bn 
that are not in our sample. 

To explore the empirical behavior of Mr-, EWN a , and 
EW Ca as a function of NIR SpT, we used a subset of 84 
M dwarfs with high S/N spectra, parallax measurements 
with precisions better than 5%, no indications of multi¬ 
plicity, and either low (EWhu <1 A, n = 57/84) or high 
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(EWhck >1 A, n = 27/84) activit y levels as evidenc ed by 
the strength of Ha emission (e.g. I West et aDl2004f ). We 
fit each of these quantities with a low-order polynomial 
as shown in Figure [TJ and defined A Mr-, AEWns, and 
AEWca as the departure from these fits at a given NIR 
SpT. For Mr- and Ca 11 , we used only the low-activity 
subset (n = 57), as t hese stars exhibit str uctural changes 
related to activity dStassun et alJ I 20 T 2 L and references 
therein) and we observe the Ca 11 feature in emission in 
some M dwarfs. We used all 84 stars for the fit to EWNa, 
because there is no clear evidence of activity measurably 
affecting the strength of the Na 1 feature in our data. 
The resulting definitions and low-order fits are: 

AEW Na = EW Na - (1-07 + 0.78 [NIR SpT] (1) 
+ 0.03 [NIR SpT] 2 ) 

AEWca = EW Ca - (5.23 - 0.93 [NIRSpT]) (2) 
AM k = M k - (4.36 -0.77 [NIRSpT]). (3) 

3. RESULTS AND DISCUSSION 

3.1. A relationship between Mk, EWN aj and [Fe/H] 

We find that a first-order linear model based on 
AEWnh and [Fe/H] yields an excellent predictor of A Mk 
(Figure 04) in the low-activity subset. The resulting 
equation is: 

A Mr- = -1.66 [Fe/H] + 0.55 AEW Na - 0.11, (4) 

which has an adjusted squared correlation coeffi¬ 
cient R 2 = 0.64 and residual scatter of 0.18 mag 
(excluding two overluminous stars that we suspect 
are binaries: 2MASS J04223199+1031188, 2MASS 
J02532611+1724324, Figure 03). This approximately 
halves the intrinsic ~ 0.35 mag scatter in Mr- from 
NIR SpT (Figured]). We explored alternat i ve definitions 
of the Na 1 features (ISchlieder et al.ir2012l : IMartin et alJ 
[Moll . and the quality of the resulting fits differs negligi¬ 
bly from that reported here. 

This model directly reveals that EWn s is sensitive to 
small changes in Mk at a given NIR SpT. This is primar¬ 
ily due to the strong sensitivity of this feature to surface 
gravity at low T e ff, a behavior that has been thoroughly 
explo red in the context of youth indicators for low-mass 
stars (ISchlieder et al.ll2012l) . With the inclusion of [Fe/H] 
in our model, we disentangle the gravity and [Fe/H] sensi¬ 
tivity of the Na 1 feature, which goes from being a binary 
indicator of youth or low gravity to having a much finer 
sensitivity to luminosity and surface gravity. Thus, this 
feature is a precise and easily-measured gravity indica¬ 
tor for M dwarfs, a stellar regime in which established 
precision grav ity measurement techniques (e.g. flicker, 
iBasticn et al.lf2C>T3h are not yet applicable. 

Our interpretation that the AMr--AEWN a-[Fe/H] re¬ 
lationship results from the gravity-sensitivity of EWN a 
depends on the observed insensitivity to gravi ty of 
our [Fe /H] measurements. Although the IMann et al.l 
( 2013a ) [Fe/H] calibrations employ gravity-sensitive fea¬ 
tures, we note that the AMR;-AEWNa-[Fe/H] relation¬ 
ship is present (with similar coefficients) with different 
indicators of [Fe/H] (e.g. the individual H and //-band 
calibrations), which are based on different features with 
varied behaviors with gravity. We also find this re¬ 
lationship using [Fe/H] derived from higher-mass com- 
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Figure 1. Top: The feature (red) and continuum (gray) definitions for the quantities used in this work. Bottom: The strengths of the 
features with NIR M subtype (a proxy for T e ff). For Na I, Ca II, and M^, we fit a low-order polynomial to these quantities for a subset 
of M dwarfs in our sample, and defined AEW]sr a , AEWc a) A as the departure from this fit (see Section n. 


Table 1 

M Dwarf Parameters and Measurements a 


2MASS ID 

[Fe/H] 

NIR SpT 

EW N a 

(A) 

EWca 

(A) 

p. b 

-‘Mnterfer. 

(Re) 

R* c 

(Re) 

i?Na d M13 Flag® 

(Re) 

00064325-0732147 

-0.11 

5.99 

7.46 ± 0.05 

0.24 + 0.07 




00085391+2050252 

0.00 

4.93 

5.38 + 0.11 

0.83 + 0.15 




00115302+2259047 

0.15 

3.62 

4.08 + 0.09 

2.38 + 0.11 




00165629+0507261 

-0.14 

4.52 

3.90 + 0.09 

1.80 + 0.11 




00182256+4401222 

-0.26 

1.89 

2.56 + 0.10 

4.45 + 0.13 

0.386 + 0.002 

0.395 ± 0.004 

0.416 + 0.027 

00182549+4401376 

-0.09 

4.10 

5.64 + 0.05 

2.22 + 0.08 




00283948-0639481 

-0.11 

4.18 

3.25 + 0.14 

2.01 + 0.15 




00313539-0552115 

-0.12 

3.68 

4.06 + 0.09 

2.12 + 0.13 




00321574+5429027 

-0.05 

4.49 

3.85 + 0.07 

1.55 + 0.10 




00383388+5127579 

-0.24 

2.83 

3.30 + 0.07 

3.25 + 0.09 





a Table ^ is published in its entirety (including parameters from the literature) 
regarding its form and content. 

Interferometric radius from fSovaiian et all <2012}) . 

C Interpolated at parallax-based M and 5 Gyr age in Dartmouth model grid. 

^ Interpolated at predicted (see Section |3. U in Dartmouth model grid. 

6 Spectrum from l^lann et al.| <2013bf) 

panions, which a r e ava ilable for 21 of our stars from 
I Valenti fc Fischer! (|2005j) . and which are immune to ef¬ 
fects within the M dwarf spectrum. We conclude that 
gravity-related systematic errors in our [Fe/H] measure¬ 
ments are insignificant. 

We considered that our NIR SpT measurements could 
similarly have systematic errors that impact our inter¬ 
pretation of the AMk- AEWNa-[Fe/H] relationship. We 
note that our NIR SpTs are based on the H 2 O-K 2 in- 


in the electronic edition of ApJ Letters, a portion is shown here for guidance 


dex, w hich i s known to correlate with T e ff to < 80K (e.g. 
iMann et all l2013bl) . By contrast, the AM^-AEWNa- 
[Fe/H] relationship extends over ~ 1.5 mag in A Mk, 
which would equate to many hundreds of degrees in T e ff. 
We conclude that, within their stated errors, our [Fe/H] 
and T e ff indicators are accurately accounting for their 
contribution to the behavior of EWNa- Therefore, the 
well-established gravity-sensitivity of EWn s is the most 
reasonable explanation for its correlation with A Mk- 
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Figure 2. (A): For 57 well-characterized nearby M dwarfs, the relationship between A M x ■ AEWns, and [Fe/H], along with iso-[Fe/H] 
contours from our model that predicts AMjf from AEWn 3 and [Fe/H], (B): For the same set of stars, the known parallax (xj-based Mx 
and our predicted Mx based on EWjj a , [Fe/H] and NIR SpT. A 1:1 relation is shown; two o verluminous targets that are likely binaries 
are highlighted. (C): For 18 M dwarfs with interferometrically-measured radii IIBovaiian ct nl .ll20 1 Hi . the fractional difference between the 
calculated (Rqsep) arid measured radii (Rbi 2 )- To derive radius, we interpolate on the Dartmouth model grid at an Mx and [Fe/H] 
appropriate to each star. We used two different values for Mx for each star: the measured parallax (7r)-based Mx (black, gray error 
bars) and the Mx predicted from our calibration of the Na I doublet (red). The R values from predicted Mx values have RMSn® ~ 9%, 
compared to RMS W ~ 5% from parallax-based Mx- We note that the mean error in Rbi 2 is ~ 2%. Both sets of radii measurements are 
consistent with < 2% offset from the Dartmouth models. GJ 725, a system with Ca II and Na I feature strengths suggestive of a-enrichment, 
is hig hlighted; the Dartm outh models with [o/Fe]= +0.2 are shown in black and red while the [a/Fe]= 0.0 results are indicated in gray and 
pink IIDottcr et al.1 (120081! explains details regarding a-enrichment in the Dartmouth models). For all other systems we used [a/Fe] = 0.0. 


This model efficiently provides Mk for M dwarfs with¬ 
out parallaxes, thereby enabling “spectroscopic distance” 
and fundamental parameter measurements of M dwarfs. 
As an illustration, we note that for one of the best- 
characterized M dwarfs in the Kepler sample, Kepler 138 
(KOI-314), we predict Mk = 5.42 ± 0.18 based on a 
single moderate-resolution NIR spectrum. This is in 
excellen t agreement with the Mk = 5.39 ± 0.25 mea¬ 
sured in lPincda et al.l (|2013li . which is based on six high- 
resolution spectra of this star and library of more than 
100 similar high-S/N spectra. Through the ti g ht rela¬ 
tions between Mr- and stellar mass (e.g. iDelfosse et al.1 
I2000D and radius, the model we present can help to 
provide determinations of M dwarf parameters and the 


location of the habitable zone for M dwarf exoplanet 
hosts. Compared wit h existing techniques for measur¬ 
ing th ese parameters dMann et al .1120 13a : iNewton et al.l 
mm, our technique benefits from a large set of cali¬ 
brators, and requires only a single, non-flux calibrated, 
moderate-resolution NIR spectrum. 

3.2. M dwarf radius from Mk 

After establishing that the Na i feature, [Fe/H], and 
NIR SpT can be used to determine M dwarf Mr-, we 
explored extending the technique to measuring M dwarf 
radius, using 18 nearby M dwarfs with interferometrically 
measured radii (Figure HP). We interpolat ed the Dart- 
mouth Stellar Evolution Program (DSEP, iDotter et al.l 
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I2008T) models at the [Fe/H] and Mk of each star to pre¬ 
dict radius, using both the parallax-based Mk and the 
spe ctroscopic Mk indicated by the techniques in Section 
13.11 We compared these r adii with the interferometric 
radii dBovaiian et al.|[2012f) . We note that the mean er¬ 
ror in the interferometrically measured radii themselves 
is approximately 2%. We found that we could derive ra¬ 
dius to ~ 5% when using the parallax-based Mk, and 
to ~ 10% with just the information available from the 
NIR spectra. We find no evidence for a > 2% systematic 
offset between Mk -based model-predicted and measured 
M dwarf radii in these 18 stars, in contr ast to interp ola¬ 
tions at LboI and T e ff (e.g. Section 5 of iBovaiian et al.1 
mm, demonstrating the advantages of using Mk and a 
possible path toward reconciling low-mass stellar models 
and observations. Using Mk to derive radius is there¬ 
fore an accurate, precise, and efficient method that can 
be applied to M dwarfs well beyond the reach of current 
interferometric facilities. 

3.3. Measurements of M dwarf a-enrichment? 

The a-elements, such as silicon and titanium, can be 
formed in Type II supernovae, and they are more abun¬ 
dant (“enriched”) relative to iron in older stellar pop¬ 
ulations. To explore the sensitivity of the Na i and 
Ca ii features to a-enrichment, we first considered a sub¬ 
set of 23 M dwarfs from our sample with a-enrichment 
([Ti/Fe] and [Si/Fe]) measur ements in th e literature from 
either direct spectroscop y dWoolf fe Wallersteinl 120051 
iChavez fe Lambert! 1200911 or based on spectrosco py of 
higher-mass companions (|Valenti &; Fischeill2005l l. We 
found that (Figure [3] A,B) a-enriched stars have higher 
(stronger absorption) AEWc a and lower (weaker absorp¬ 
tion) AEWNaj as expected since calcium is an a-element 
and sodium is not. Moreover, for the subset of our M 
dwarfs where measurements are available, high AEWc a 
and low AEWn s are also consistent with kinematic and 
color-based indicators of old age. A Toomre diagram 
(Figure [3] C, ISandage fe Foutslll987fl . which shows stellar 
Galactocentric vertical (U 2 ), radial (W 2 ), and rotational 
(V 2 ) kinetic energies, can dynamically separate older and 
younger stars: older stars have a wide distribution of ki¬ 
netic energies, while younger stars have lower kinetic en¬ 
ergies. The lower metallicities of older stars also separate 
these stars well in the (J — H, H — K) color-color plane 
(Figure EH D) , due to H~ continuum and H 2 0 opacity 
effects f Leggettl 11992 ). These observational indications 
of abundance and age suggest that the Na i and Ca ii 
features jointly indicate a-enrichment in M dwarfs. 

Figure [4] shows the plane defined by AEWn 3 and 
AEWca- We note that, for the lowest AEWn 3 val¬ 
ues, targets with low and high Ha are generally well- 
separated by AEW Ca (with the exception of CM Dra, an 
active close binary system). This is consistent with the 
behavior of AEWca being dominated by chromospheric 
activity, which is supported by our observati o n of Ca ii in 
emission and the findings of lKafka fc Honevcuttl (2006). 
We also note tha t CM Dra, a system that is suspected 
of a-enrichment (|Feiden fc Chaboverll2014!j h stands out 
in this plane, as does the binary system GJ 725. In both 
systems, models with [a/Fe]=0.2 provide a better match 
to observations than models with [a/Fe]=0.0 (e.g. Fig- 
ure[2|). We posit a simple interpretation of these features: 


stars with low AEWn 3 are likely either low-gravity or a- 
enriched, and the activity-sensitivity of AEWca can dif¬ 
ferentiate between these younger (more active) or older 
(less active) populations. As we only probe a limited 
range of [a/Fe] and [Fe/H], this interpretation is not con¬ 
clusive, but is strongly suggestive and provides a promis¬ 
ing direction for further exploration and the development 
of better techniques for measuring M dwarf compositions. 

4. CONCLUSION 

The Ca n and Na i spectroscopic indicators provide 
a wealth of new and easily obtained information about 
M dwarf stellar properties. The broad applicability 
of these physically-motivated indicators and associated 
techniques in mid- to late-type M dwarfs, their possi¬ 
ble constraints on a-enrichment, and their simplicity of¬ 
fer a substantial improvement over existing methods for 
measuring M dwarf parameters. These indicators and 
techniques will be useful for appraising M dwarf plane¬ 
tary systems discovered in the coming years (by e.g. by 
TESS and K2 ), and for prioritizing nearby M dwarf ex¬ 
oplanet systems for planetary atmospheric characteriza¬ 
tion with JWST. By providing stellar mass, luminosity, 
and radius through Mk, they will provide a necessary 
foothold for measuring planetary parameters and the lo¬ 
cation of the habitable zone. They also provide a path 
toward more precise distances to Galactic M dwarfs, for 
which the current limiting factor is the scatter in the 
color-luminosity relati ons used to estimate distance (e.g. 
iBochanski et al.l[2013l) . The Sloan Digital Sky Survey 
bandpass contains the Ca n and Na i /-ba nd features, 
as well as indicators of m etallicity and T e ff dWest et al.l 
1201 It iDhital et al.l 1201211 , potentially enabling archival 
analysis with the Na i and Ca n indicators. Finally, we 
note that the Ca n feature will be observed by the Gaia 
RV Spectrometer and may provide leverage for the mea¬ 
surement of a-enrichment and activity levels in the large 
sample of M dwarfs observed by Gaia. The measurement 
of the features described here is relatively observationally 
inexpensive: a single R ~ 2000 NIR spectrum is suffi¬ 
cient, and is substantially easier to obtain than the paral¬ 
lax measurements, high-resolution spectra, or broadband 
flux-calibrated spectra required by other techniques. 
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Figure 3. The behavior of the Ca II and Na I features with a-enrichment. (A,B): For 23 M dwarfs, the reported [a/Fe] parameter is 
the weighted averages of literature measurements of [Si/Fe] and [Ti/Fe] of a higher-mass common proper motion companion or the M 
dwarf itself, as a function of AEWNa and AEWca- High AEWca and low AEWNa are clearly associated with a-enrichment. (C): The 
Toomre diagram for M dwarfs with well-measured UVW velocities, and nominal Galactic population boundaries. Stars that are dynamicall y 
consistent with being old also have high AEWca- (D): The color-color diagram with delineations of different populations (|Leggetti119921 b 
M dwarfs with high AEWca have colors consistent with old ages. 
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Figure 4. The AEWca and AEWNa values for M dwarfs in our sample, color-coded by their level of activity as indicated by Ha. We 
posit a simple interpretation: stars with low AEW^a are likely either low-gravity and young, or a-enriched and old, and AEWc a can 
discriminate b etween these two due to fi lling-in related to chromospheric activity Two M d warf systems (CM D ra and GJ 725) with possible 
a-enrichment JFeiden &; ChaboveiH201^) are highlighted. Measurements on spectra from JMann et al.ll2013bl) are indicated separately and 
have been corrected for the use of a different instrument (0.96 A offset). Mean errors for the stars observed in our program are indicated 
by the error bars in the lower left. 
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